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' Abstract. We present the analysis of the observations of standard stars for 

linear polarization obtained from 1999 to 2005 within the context of the cali- 
bration plan of the FORS1 instrument of the ESO VLT. We have considered 
observations carried out both in imaging polarimetric and in spectropolarimet- 
ric mode. Broadband polarization was obtained in the Bessel BVRI filters; 
spectropolarimctry was obtained with various grisms covering different optical 
\q \ ranges and with a typical resolution of a few hundreds. Spectropolarimetric data 

. have been convolved with the transmission functions of the Bessel filters, which 

enabled us to calculate polarization values equivalent to broadband polarization 
\Q ] measurements in imaging mode. Finally, for each star, instrument mode, and 

for each Bessel filter band, we have calculated an average polarization value 
obtained after filtering all available data with a fcc-clipping algorithm. 
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When taking polarimetric measurements of scientific targets, observations of 
standard stars for linear polarization are of crucial importance to calibrate and 
monitor instru ment performances. Popular lists o f standard stars are, e.g., those 
by Serkowski (|l974l ) and Hsu and Breger (|l99£l ). Unfortunately, most of the 
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'classical' standard stars for polarization presented in these works are too bright 
for the instruments of the large size class telescopes. 

At the ESO Very Large Telescope (VLT), observations of standard stars 
deemed to exhibit either zero polarization or substantial polarization are rou- 
tinely performed within the context of the FORS1 calibration plan. Usually, 
a star known to exhibit a large signal of linear polarization is observed during 
those nights when a science target is observed in polarimetric mode. Occasion- 
ally, a non polarized star is also observed. However, the stability of these stars 
used as standard for linear polarization has never been extensively checked so 
far. Therefore, we have decided to retrieve from the archive a large sample of 
observations of polarimetric standard stars obtained with FORS1 and to check 
their consistency. 
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2. Instrument and observations 

F0RS1 is the visual and near ultraviolet FOcal Reducer and low dispersion 
Spectrograph mounted at the Cassegrain focus of one of the four 8 m units of 
the ESO VLT, and works in the wavelength range 330 — 1100 nm. FORS1 is 
equipped with polarization analyzing optics, which include a half-wave and a 
quarter-wave retarder plates, both superachromatic and a Wollaston prism with 
a beam divergence of 22". This allows the measurement of linear and circular 
polarization, both in direct imaging and spectroscopy. In imaging polarimetric 
mode (IPOL) the field of view is 6.8' x 6.8' and the magnitude limit is R = 23, 
for a 1 % accuracy in the polarization measure, and with 1 h exposure time. In 
spectropolarimetric mode (PMOS), the magnitude limit is between R = 17.2 
and 19.3, for a spectral resolution between 260 and 1700 (with a 1" slit width, 
depending on the grism inserted). The instrument is described in Appenzeller 

et al. (1998). 

We have retrieved from the ESO data archive (http://archive.eso.org) 



all observations of polarimetric standard stars taken with FORS1 from April 
1999 to March 2005. A preliminary inspection to the data was performed to 
identify and discard a few saturated exposures and some frames with poor im- 
age quality. The list of the useful observations is presented in Tabled About 210 
observations in IPOL mode were obtained with the BVRI Bessel filters (obser- 
vations with the U Bessel filter are not possible in polarimetric mode as the filter 
is situated in the same wheel as the Wollaston prism). About 130 observations 
in PMOS mode were obtained with a variety of grism and filter combinations. 
Both for IPOL and PMOS modes, the linear polarization was usually measured 
taking a series of four observations with the half wave retarder waveplate at 
a = 0°, 22.5°, 45°, and 67.5°, respectively, where a indicates the angle between 
the acceptance axis of the ordinary beam of the Wollaston prism and the fast 
axis of the retarder waveplate. The acceptance axis of the Wollaston prism was 
always aligned to the North Celestial Meridian, which therefore represents the 
reference direction for all linear polarization measurements of this paper. 



3. Data reduction 



Stokes Q and U parameters are defined as in Landi Degl'Innocenti et al. (2007) 



with the reference axis coinciding with the North Celestial Meridian. In the 
following, we will consider the ratios Q/I and U/I, adopting the notation 



Pq = | and Pu = j. (1) 



Pq and P\j were measured by combining the photon counts (background sub- 
tracted) of ordinary and extra-ordinary beams (/° and / e , respectively) observed 
at retarder waveplate positions a = 0°, 22.5°, 45°, and 67.5°, as given by the 
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following formula: 



Pq 
Pu 



i J ( r-r \ _ ( f°-.r \ 
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a=45° 



/°-/ c 
/°+/ c 



a=22.5 c 



_ f f°-f e \ 
\f°+f c )a-- 



(2) 



:67.5° 



(see FORS1/2 User manual, VLT-MAN-ESO-13100-1543). The error on P Q or 
Pu is 

2 „ / .„ x2 



a /° + V(/°+/°)^ / c 



((/° 



+f e y 



f° 
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a=</>0 



a=45°+0o 



(3) 



where 4>q = 0° if X = Q and </>o = 22.5° if X = U. If the polarization of the 
targets is always small, in order to give an estimate of the quantities op x one 
can substitute in Eq. © f e = f° = f, with / independent of a. Also, if one 
can assume ajc = a jo = ot we have 



°7 

°Px = -T ■ 



(4) 



If the error on the photon counts is entirely due to photon-counts, we get 

1 1 



(5) 



From Pq and Pjj we have obtained the total fraction of linear polarization 
Pl and the position angle 9 (see Landi Degl'Innocenti et al. 2007). For the 
cases where polarimetric errors are small with respect to the signal (cp Q <C Pl, 
ap v <C Pl), the errors on Pl and 9 are 



cos 2 (2(9)cj| 3 + sm 2 (29)a 2 P 



1/2 



(6) 



On = - 

2 



Note that if ap Q = op v one gets 



sm 2 (29)a 2 PQ + cos 2 {29)a 2 Pu 



1/2 



Pl 



and 



ap L = u Pq = a Pu 



1 ap L 



2 Pl 



(7) 

(8) 
(9) 
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3.1. IPOL data 

All the science frames were bias subtracted using the corresponding master bias 
obtained from a series of five frames taken the morning after the observations. No 
flat-field correction was carried out as it is irrelevant for the purpose of measuring 
the polarization through Eq. ©. The flux in the ordinary and extra-ordinary 
beams was measured via simple aperture photometry, that was performed using 
the DAOPHDT package implemented in IRAF. Once the fluxes for the ordinary 
and extraordinary beams for each retarder waveplate position were measured, 
we used a dedicated C routine to calculate P'q and P{j through Eq. ©, the 
corresponding errors via Eq. (J3J), then the fraction of linear polarization P^ 
and th e position angle 8' with Eqs. (6) and (7) of Landi Degl'Innocenti et al. 
(2007), with the corresponding errors given in Eqs. © and ((7J) of this paper. To 
compensate a chromatism problem of the half wave plate (see Sect. 3.2) a new 
position angle 9 was obtained as 

= 0'- ee (F) (10) 

where e#(F) is a correction factor that depends upon the filter F. The e#(F) 
values are tabulated in the FORS1/2 user manual (see also 
http://www.eso.org/instruments/fors/inst/pola.html). We obtained the 
final Pq and P\j values as 

P Q = Plcos{29) 

P v = Pi sin(2#) . [Ll) 
with 6 given by Eq. (note that P L = P[). 

3.2. PMOS data 

PMOS data were pre-processed with the packages for spectra analysis imple- 
mented in IRAF. Spectra in the ordinary and extraordinary beams were bias- 
subtracted, then optimally extracted and wavelength calibrated using IRAF 
routines, and finally processed using dedicated C routines to calculate Pq, Pu, 
Pl, and 6 with the corresponding errors using Eqs. © and (j3J) of this paper, 
and Eqs. (6) and (7) of Landi Degl'Innocenti et al. (2007]). 

Figure Q shows the results for the star Ve 6-23 observed in the night 17 
December 2003 with the 1501 grism and the GG435 filter. The position angle 
of interstellar polarization is expected to have a constant value independent of 
wavelength. The botton right panel of Fig. ^ shows that this is not the case for 
the star Ve 6-23 (as well as for all other stars of our sample). This is due to 
a chromatism problem of the half-wave retarder waveplate already discussed in 
Sect. 3.1. 

To compare the IPOL with the PMOS observations we convolved the po- 
larized spectra obtained in PMOS with the transmission functions of the BVRI 
Bessel filters used in IPOL mode. For each filter F we calculated 

f~d\P Q (\)I Q (\)T F (\) _ / oo dAP^(A)/ c/ (A)T F (A) 

rQ[ "> / °°dA/Q(A)T F (A) / °°dA%(A)T F (A) [ ^ 
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6000 7000 8000 9000 rt 6000 7000 8000 9000 



Wavelength (A) Wavelength (A) 

Figure 1. Star Ve 6-23: Pq, P' v , Pl,, and 9' observed in the night 
17 December 2003 with grism 1501 and filter GG435. The non- 
constant 9 value is due to the chromatism of the half waveplate. 
The thin solid line shows the offset angle eg (A) (data are available at 
http://www.eso.org/instruments/fors/inst/pola.html) that has to be 
subtracted to the observed position angle to compensate for the waveplate 
chromatism (in the figure, eg (A) is offset by a constant to allows its visualiza- 
tion) . 

where Tf is the transmission function of the F filter, and 

iq = cr + n u=o° + (/° + n u=45° 

IU = (f° + n U=22.5° + (/° + / G ) U=67.5° • 

The polarization values so obtained have been finally modified according to the 
procedure followed for the IPOL data, in order to compensate for the chromatism 
of the half wave plate. 

Error bars obtained in PMOS mode are smaller than those obtained in 
IPOL mode. This is a consequence of the fact that the signal-to-noise ratio of 
the observations is basically limited by the full well capacity of the CCD pix- 
els, or by the hardware limitations of the digital-analogic converter. Rebinning 
spectropolarimetric data permits one to integrate the signal over a much higher 
number of pixels than possible in imaging mode, leading to a much higher signal- 
to-noise ratio. A comparison between IPOL and rebinned PMOS data can be 
seen in Figs. |U and for stars Ve 6-23 and CD —28 13479, respectively, in the B 
band. 
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Figure 2. IPOL data (full circles) and rebinned PMOS data (open circles) 
for the star Ve 6-23 in the B band 



3 § 



5 - 



1500 



2000 2500 
JD-2450000 



IS -J 

- 1 

D 



1500 2000 2500 3000 
JD-2450000 



3 § 



5 - 



1500 2000 2500 3000 
JD-2450000 



t 5 

C 180 - 1 



1500 2000 2500 3000 
JD-2450000 



Figure 3. IPOL data (full circles) and rebinned PMOS data (open circles) 
for the star CD -28 13479 in the B band 
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4. Results and discussion 

Figures [2 and |3] show the polarization observed in IPOL and in PMOS mode 
for stars Ve 6-23 and CD —28 13479, respectively, plotted as a function of the 
observation epoch. The plot scales are the same for both stars. The data for the 
star Ve 6-23 appear much more scattered than those for the star CD —28 13479. 
Whereas some scattering may well be due to undetected instrument or data 
reduction problems, Fig. 21 suggests that the polarization of the star Ve 6-23 
may exhibit a short-term variability. 

For each star, we have grouped all observations obtained with similar in- 
strument mode and with the same filter (or convolved with the same trasmission 
function, in case of PMOS observations). From each group so obtained we have 
calculated the medians of the observed Pq and Pjj values, Px (with X = Q 
and U), and the median absolute deviations (MAD), i.e., the medians of the 
distributions 

\{Px)i-Px\. (14) 

Setting o = 1.48 MAD (e.g., Huber 1981, pp. 107-108), we have then rejected 
those (Px)i values for which 

\{Px)i-Px\ >3a (15) 

Finally, from the remaining values, we have calculated the weighted averages Pq 
and Pjj: 

Px = • (16) 

To each average value, we have associated the error given by 



({Px)i 

E, 



2 

-Px' 



The results are shown in Tabled (for the stars with high polarization signal) 
and |31 (for the stars with low polarization signal). Stars observed only once were 
not included in the tables. For stars with less than four observations we did not 
run the fee-clipping algorithm. 

4.1. Polarized stars 

Table [2] may be used as a reference to check instrument performance and sta- 
bility within the quoted error bars, both in IPOL and PMOS mode. Note that 
relatively large errors may point to a star's instrisic variability, especially if large 
errors are associated to large data sets. Observations within similar filter bands 
are fairly, but not always fully consistent, see Sect. 4.2 below. 

It should be noted that the actual correction that has to be applied to the 
broadband polarization measurements to compensate for the chromatism of the 
retarder waveplate depends on the shape of the star's spectral energy distribution 



Standard stars for linear polarization 



9 



(convolved with the transmission of the telescope optics). In fact, we have 
applied a correction that is independent of the star's colour. This is probably 
the reason why the position angles of the stars of Table [21 are slightly filter 
dependent. For the same reason, caution should be adopted when comparing 
the results reported in Table [21 with polarimetric observations obtained with 
other instruments. 

4.2. Unpolarized stars 

All stars observed in IPOL mode and reported in Table [31 have Pq and Pjj 
values consistent with zero. This means that all these stars may be considered as 
unpolarized standard stars within a typical accuracy better than 3 x 10 -4 . At the 
same time, the available IPOL observations do not show evidence for significant 
instrumental polarization in the center of the instrument field of view 1 . 

WD 2007-303, that was observed in PMOS mode only, is polarized at about 
0.5% level in Pq, and should not be used as unpolarized standard star. 

The remaining PMOS observations hint that there is a small instrumental 
offset in Pq . The weigthed average of all Pq values of Table [3] (but without 
considering the star WD 2007-303) is Pq(B)\ = 0.07 ± 0.01 %, and Pq{V)\q = 
0.09 db 0.01 % for the B and the V filter, respectively (the errors were calculated 
with Eq. (|17|> ). The averages of all Py values of Table is fully consistent with 
(with an error bar of about 0.005 %) both in the B and in the V filters. The 
observed Pq offset, that appears instrinsic to the PMOS mode only, may be 
associated to some but not all grism+filter combinations, and deserves further 
investigation. 

Acknowledgments. L. Fossati acknowledges ESO DGDF for a four month 
studentship at ESO Santiago/Vitacura. 
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Table 2.: Weighted averages of the observations of polarized stars 
in IPOL and PMOS mode with various instrument settings. The 
averages were calculated after running the fctr-clipping procedure 
explained in the text. N s and N c represent the number of data 
before and after running the fcc-clipping procedure, respectively. 
Note that for N < 3 no fccr-clipping procedure was applied. 
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IPOL 


v 


4.85±0.02 


0.41±0.07 


4.87±0.02 


2.40±0.38 


3 


3 




IPOL 


R 


4.69±0.03 


0.55±0.02 


4.72±0.03 


3.36±0.12 


3 


3 




IPOL 


I 


4.32±0.01 


0.40±0.04 


4.34±0.01 


2.65±0.26 


2 


2 




PMOS 


B 


4.62±0.02 


0.57±0.01 


4.66±0.02 


3.54±0.06 


4 


4 




PMOS 


V 


4.97±0.01 


0.51±0.02 


5.00±0.01 


2.94±0.10 


3 


3 


BD-14 4922 


IPOL 


B 


-0.98±0.03 


5.63±0.01 


5.71±0.01 


49.93±0.15 


4 


4 




IPOL 


V 


-0.99±0.02 


6.07±0.02 


6.15±0.02 


49.62±0.08 


5 


5 




IPOL 


R 


-0.98±0.05 


5.69±0.07 


5.77±0.07 


49.90±0.24 


2 


2 




PMOS 


B 


-0.97±0.01 


5.66±0.01 


5.74±0.01 


49.88±0.07 


5 


4 




PMOS 


V 


-0.93±0.01 


6.06±0.01 


6.13±0.01 


49.35±0.03 


5 


3 




PMOS 


I 


-0.70±0.03 


5.07±0.01 


5.12±0.01 


48.92±0.19 


2 


2 


BD-12 5133 


IPOL 


B 


1.87±0.04 


-3.95±0.05 


4.37±0.05 


147.64±0.28 


6 


5 




IPOL 


V 


1.75±0.04 


-4.00±0.04 


4.37±0.04 


146.84±0.25 


5 


4 




IPOL 


R 


1.63±0.02 


-3.68±0.02 


4.02±0.02 


146.97±0.13 


16 


13 




IPOL 


I 


1.10±0.29 


-3.39±0.07 


3.57±0.09 


143.99±2.27 


3 


3 




PMOS 


B 


1.94±0.04 


-3.94±0.03 


4.39±0.03 


148.09±0.23 


5 


4 




PMOS 


V 


1.80±0.04 


-4.05±0.01 


4.43±0.02 


146.99±0.22 


6 


4 




PMOS 


I 


1.03±0.47 


-3.37±0.06 


3.53±0.11 


143.45±3.73 


2 


2 


BD-13 5073 


IPOL 


R 


2.10±0.01 


-2.99±0.03 


3.66±0.02 


152.55±0.11 


2 


2 




PMOS 


V 


3.40±1.59 


-2.64±1.17 


4.30±1.47 


161.11±8.77 


2 


2 



Standard stars for linear polarization 



Table 3.: Weighted averages of the observations of unpolarized 
stars in IPOL and PMOS mode with various instrument settings. 
The meaning of the various columns is the same as in Table |3 



Star 


KTrn-lr- 




RflnH Pn (%1 
JDailLl JTQ \ /0 ) 


Pu (%\ N 


N 




WD 2359-434 


IPOL 


R 


0.03±0.06 


-0.20±0.09 


4 


4 




IPOL 


1 


-0.00±0.01 


0.03±0.04 


4 


4 


HD 10038 


PMOS 


13 


0.07±0.01 


-0.03±0.02 


3 


3 




PMOS 


V 


O.lliO.Ol 


-0.03±0.01 


2 


2 


HD 13588 


PMOS 


B 


0.10±0.01 


0.00±0.01 


4 


4 




PMOS 


V 


0.10±0.01 


0.00±0.01 


4 


4 


VV VJ UOIU-DOO 


TDflT 
lr ULi 


V 


n ni -i-n c\a 

U.UllCU.U'rt 


— U.UOlCU. 1U 


Q 
D 


Q 
U 




IPOL 


R 


0.04±0.09 


0.01±0.14 


2 


2 


HD 42078 


PMOS 


B 


0.04±0.01 


-0.01i0.01 


4 


2 




PMOS 


V 


0.07d=0.01 


-0.02±0.01 


6 


6 


HD 64299 


PMOS 


B 


-0.02±0.12 


-0.04±0.04 


2 


2 




PMOS 


V 


-0.05±0.13 


-0.03±0.04 


2 


2 


WD 0752-676 


IPOL 


13 


-0.00i0.01 


0.06±0.03 


4 


3 




IPOL 


V 


0.00±0.01 


0.03±0.05 


3 


3 


HD 97689 


PMOS 


B 


0.16±0.04 


O.OliO.Ol 


3 


3 




PMOS 


V 


0.14±0.08 


-0.00i0.01 


2 


2 


WD 1615-154 


IPOL 


V 


0.03±0.24 


-0.05±0.03 


2 


2 




IPOL 


R 


-0.01±0.04 


-0.02±0.12 


2 


2 




PMOS 


13 


0.07±0.02 


-0.03±0.02 


5 


5 




PMOS 


V 


0.13±0.01 


-0.02±0.02 


4 


4 


WD 1620-391 


IPOL 


B 


0.02±0.03 


-0.02±0.02 


3 


3 




IPOL 


V 


-0.01iO.01 


-0.01i0.01 


1 


3 




IPOL 


R 


-0.00±0.01 


0.01±0.01 


5 


4 




IPOL 


1 


0.01±0.01 


O.OliO.Ol 


4 


2 




PMOS 


13 


0.06±0.02 


0.00±0.01 


2 


2 




PMOS 


V 


0.08±0.02 


0.00±0.01 


2 


2 


HD 176425 


PMOS 


V 


0.15±0.13 


-0.08±0.06 


2 


2 


WD 2007-303 


PMOS 


B 


0.44±0.01 


-0.09±0.01 


4 


3 




PMOS 


V 


0.51±0.02 


-0.05±0.01 


4 


4 


WD 2039-202 


IPOL 


R 


0.19±0.23 


-0.14±0.40 


2 


2 




PMOS 


13 


0.03±0.03 


0.00±0.01 


2 


2 




PMOS 


V 


0.03±0.01 


-0.00i0.01 


2 


2 




PMOS 


1 


0.10±0.01 


O.OliO.Ol 


2 


2 


WD 2149+021 


IPOL 


V 


0.01±0.03 


-0.04±0.01 


2 


2 




IPOL 


R 


0.14±0.14 


-0.02±0.09 


3 


3 




PMOS 


V 


0.09±0.01 


-0.01±0.03 


2 


2 




PMOS 


1 


0.16±0.01 


-0.01±0.04 


2 


2 



